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After recalling some connections between the Spontaneous Breakdown of Chiral 
Symmetry (SB^S) and the spectrum of the Dirac operator for Euclidean QCD on a 
torus, we use this tool to reconsider two related issues : the Zweig rule violation in 
the scalar channel and Vine dependence of SB%S order parameters on the number 
Nt of massless flavoursEI. The latter would result into a great variety of SBxS 
patterns in the (Nf,N c ) plane, which could be studied through Leutwyleg-Smilga 
sum rules in association with lattice computations of the Dirac spectrum!]. 



1 Chiral symmetry and its breakdown 

At low energies, QCD cannot be described in a perturbative way : its 
degrees of freedom (quarks and gluons) are different from its asymptotic 
states (hadrons). To analyse the theory, one can fortunately follow an- 
other path and consider its symmetries. If we set to zero the masses of 
the Nf lightest quarks (chiral limit), the right- and left-handed fermions can 
be rotated separately in flavour space and QCD exhibits a chiral symmetry 
SUL(Nf) x SUfi(Nf) x LV(1). This symmetry is spontaneously broken along 
its axial directions by the fundamental state of the theory, generating pseu- 
doscalar Goldstone bosons. When light quark masses are turned on, they 
appear in the QCD Lagrangian as small perturbations breaking down explic- 
itly chiral symmetry and providing the Goldstone bosons a light mass. 

In Nature, Nf — 3 flavours are light enough in comparison with QCD 
scale (~ 1 GeV) for this treatment, and the Goldstone bosons are identified 
with the light pseudoscalar octet (ir,K,r)). It is often stated that SBxS is 
triggered by a large condensation of quark-antiquark pairs in the vacuum: 
(0|gg|0) = — Y,(Nf), but this assumption has to.(and is about to) be tested 
experimentally, at least in the two-flavour caseu. Supposing the leadership 
of the quark condensate in SBxS, the Standard Chiral Perturbation The- 
ory (SxPT)o provides an effective description of the pseudoscalar-octet phe- 
nomenology. With the same goal, the Generalized xPTO keeps the quark 
condensate as a free (and possibly small) parameter and, at a fixed order of 
accuracy, takes into account contributions neglected by SxPT. 
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We highlight the particular status of the two-point correlator 

E^(q 2 )S^ ^ijd+x e***(Sl\T{V*(x)Vj(0) - 4(x)4(0)}|O) (1) 

where V (A) are the Noether currents for vector (axial) rotations in flavour 
space, and is the vacuum of the massive theory. E^„ is a chiral order 
parameter for any momentum q. In the chiral limit, a massless pole arises 
because of the Goldstone bosons coupling to the axial current. £ M „(0) tends 
then to —rj fJi ^F 2 (Nf), yielding the decay constant of the Goldstone bosons 
for Nf massless flavours. A non- vanishing 2^(0) is thus a necessary and 
sufficient condition for SB^S, by contrast with other order parameters, whose 
vanishing is not a clear signal for the full symmetry restoration. 

2 iV/-dependence of SB^S order parameters 

A common approximation leads to neglect the (supposed) weak in- 
dependence of SBxS order parameters due to light quark loops: this effect is 
suppressed in the large- N c limit_and violates the Zweig rule, both considered 
as good approximations of QCDEl. On the other hand, the scalar channel does 
not obey large- iV c predictions, and all chiral order parameters should vanish 
for Nf/N c large enough, where the perturbative QCD /3-function predicts the 
end of confinement. Furthermore, recent lattice simulations!] find large varia- 
tions in the SB^S pattern when Nf increases, and a sum rule estimated, based 
on experimental data in the scalar channel, asserts H(Nf) could decrease by 
half from 2 to 3 flavours. 

To study this problem, Euclidean QCD on a torus (of volume V = L 4 ) 
seems an appealing tool, because SB%S is related to the lower end of the 
spectrum of the Dirac operator: H[G] = 7^(9^ + iG^). For each gauge 
configuration, this Hermitean operator admits only real eigenvalues (e.v.s) 
A„. The spectrum is symmetrical with respect to zero, because {H, 75} = 0, 
and the winding number u[G] of the gauge configuration provides the number 
of exactly zero e.v.s. We label positive e.v.s in the ascending order with a 
positive integer, and denote the negative part of the spectrum A_„ = — A n . 

In a d-dimejisional space, Dirac e.v.s turn out to be uniformly bound by 
the free theory!!: |A„[G]| < u) n = C ■ n x l d /L, where C depends only on the 
geometry of the space-time manifold, but neither on G, n nor V = L d . This 
bound shows the paramagnetic response of the Dirac spectrum to an external 
gauge field. We see the lowest e.v.s accumulate around zero for L — > 00. 

The eigenvalues A and the eigenvectors 4>(x) of the Dirac operator are 
independent of quark flavours and quark masses, and they can be used to 
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project and integrate fermionic fields in any functional integral. The following 
propagators and the (normalized) determinants are involved: 

n J L J n>0 

In particular, the fermionic determinant reads: det(A7 — iH) = ■ A(rrij\G). 
In the chiral limit toi = ._. = mjv, = to ^ (the other quark masses 
remaining fixed), (qq) yieldscl: 

E(iV» =lim^4 < dx Tr S , 1 (x,x|G)> W/ = hm— <C^ ^2 - A2 > jV/ 

n n 

(3) 

where lim means the limit L — > oo, followed by to — > 0. The average over 
gauge configurations is: 

« r >at /= Z- 1 f dn[G\ e- s[G] r A^(to|G) J[ A(to j |G), (4) 

j>N f 

where the dependence on Nf is explicit. This average is conveniently nor- 
malized (<C 1 ^>N f = 1) Emd requires a non-perturbative mgularization and 
renormalization. In a similar way, (^) once projected yieldalB 

o , 1 v— * to to 

with the transition probability: J^ n — \ Jdx 4>\{x)^ tIl (j) n {x)\ 2 . In the 

chiral limit m —> 0, F 2 and S are essentially dominated by the lowest Dirac 
e.v.s, because of the factors to/(to 2 + A 2 ). The quark condensate is only 
sensitive to the e.v.s accumulating like 1/i 4 in the gluonic average, while 
1/L 2 is a sufficient speed for an e.v. to contribute to F 2 . 

In A(to|G), infrared and ultraviolet e.v.s can be split with respect to a 
cutoff A. If K is the integer such as wr — A, we have 

A = toI"IA ir (to|G)A uv (to|G), A IR (TO|G) = n m2 2 +A2fc[ 2 G] < 1. (6) 

1 A TO + LOT 

k=l k 

Ayv needs regularization and renormalization. However, in (^), order pa- 
rameters dominated by the lowest Dirac e.v.s should be essentially sensitive 
to Air, which is bounded by 1 and more and more efficiently suppressed 
when Nf goes up. Infrared-dominated order parameters should then decrease: 
E(Nf + 1) < S(iV/) and F 2 (N f + 1) < F 2 (N f ). Moreover, since the quark 
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condensate is more sensitive to the lowest e.v.s than F 2 , its suppression should 
be stronger. 

The dependence of £(iV/) on the (Nf + l)-th quark, called s, is given by 
the connected part of the correlator: 

- — E(JV» = lim / dx(0\uu(x) ss(0)|0) c = Tl z (m s ) (7) 
om s m— >0 J 

so that T,(Nf) = T,(Nf + 1) + f™ s dfi Hz(^)- If the number of flavours lies 
just below n cl -it(N c ), the large difference E(JV/) — S(iV/ + 1) implies, through 
Hz, a significant Zweig rule violation in the ++ channel which is actually 
observed for Nf = 2 — 3. Besides, the large- N c limit throws n cr jt to infinity: 
if «crit(3) lies near 3, the 1/N C expansion should converge very slowly to the 
phenomenology in the scalar sector. The study of chiral phase transitions, far 
from being academic, could thus deepen our understanding of actual QCD. 



3 Leutwyler-Smilga sum rules : extensions and applications 

In the (Nf,N c ) plane, SB^S order parameters could therefore exhibit very 
different behaviours for their decrease and their vanishing, leading to a rich 
chiral phase structure for QCD. The lowest Dirac eigenvalues may shed some 
light col. this problem, through sum rules first described by Leutwyler and 
SmilgaEa. In a large torus, the lightest excitations (the pseudoscalar mesons) 
dominate the partition function of Euclidean QCD, because heavier states are 
exponentially suppressed. Moreover, since periodic boundary conditions are 
chosen, the effective Lagrangian density £ e ff describing the Golstone bosons 
in a box with periodic boundary conditions is identical to its equivalent in 
an infinite volumai3 : neither Lorentz-breaking terms, nor volume-dependent 
coefficients arise, once on the torus. 

The QCD and xPT representations of the partition function can thus be 
matched for large volumes 

djjL[G] e~ s[G] det(-ip+ M) - / dfi[U] exp \-S cS (U, dU, Me ie > N f) \ , (8) 



where U(x) G SU(iV/) collects the pseudo-Goldstone bosons, M = (1 — 
75)Af/2 + (I + 7s)M'l'/2 is the light quark mass matrix (with positive real 
eigenvalues), and the vacuum angle, conjugate to the winding number v[G] 
in the QCD Lagrangian. (||) leads to sum rules for {(Jk)v'- 

^ = Ei ^ = A » I d ^ e ~ SlG] ( II X2 n) Nf ^ ( 9 ) 
n>0 « Jv «>0 
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where the functional integral is performed over the gluonic configurations 
with a fixed winding number v, and A v is a normalization constant such that 

<1>« = 1- 

The sums {&k)v are especially sensitive to the lowest Dirac eigenvalues, 
and should therefore reflect SB^S. Actually, Leutwyler and Smilga obtained 
the asymptotic behaviour of such sums for L — > oo, for instancdl3 

{ 2h m + w\y { )v i6[Nf+Mm+w\r-iy (1U) 

Two questions arise: how would these limits change in a phase where chiral 
symmetry is broken but (qq) vanishes ? what are the finite-size corrections to 
these asymptotic results ? 

( |l0| ) were derived by considering only the leading order of S e g in S^PT. 
The finite-size corrections to the sum rules are due to higher orders and should 
be hard to disentangle even at intermediate volumes, if their contributions are 
small compared to the one due to the quark condensate. On the other hand, 
a small (aq) leave space for large finite-size effects. G%PT is designed for this 
situation^ and e.g. for 04 at v = 0, one ends up with the sum rulea: 

(ff4) ° = i 6jV/ (^~l) { (yS]2)2 + [VY?) ' 4i?2[3Zs _ Zp ~ N f A] 

+4F 4 [3(Z 2 S + Z 2 P ) - 2Z S Z P + A 2 - 2N f (Z s + Z P )A]} + . . . (fl) 

where A, Z$ and Zp are chiral order parameters arising in the effective La- 
grangian. The dots remind of higher order terms: ( |iT|) neglects contributions 
that are supposed to remain small for large volumes. If S = 0, the asymp- 
totic behaviour of (tJ4)o dramatically changes from V 4 to V 2 , indicating a 
phase transitional. On the contrary, if the quark condensate is small but 
non- vanishing, ( |l l| ) leads, even at intermediate volumes, to a large departure 
from (|l^) , and its shape is related to the detailed pattern of chiral symmetry 
breaking (F 2 , S, A, Z s ...). 

These sum rules might be exploited through lattice simulations, since : 

1. the A^y-dependence is explicitly displayed in the gluonic average (||), 

2. renormalization group-invariant ratios of the sums Uk can be studied, 

3. computing the sums requires a set of gluonic configurations (v = 0) 
and the resulting lowest Dirac e.v.s, obtained by diagonalizing the discrete 
version of Tp — D 2 + to avoid the problematic doublers from Jj). 



a This result agrees with the discussion in Sec. ^ : the l/L 2 -eigenvalues contribute to F 2 
(SB^S), but not to S. Such a situation would take place for Nt > n cr ; t (Af c ). 
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Once these sums are computed for various (large) volumes, fits with sum 
rules similar to ( pi] ) could determine chiral order parameters for the desired 
point in the (Nf, N c ) plane, until a full restoration of chiral symmetry. 

4 Conclusion 

Euclidean QCD is a powerful tool to investigate chiral symmetry, since its 
spontaneous breaking is reflected by the lowest eigenvalues of the Dirac 
operator. This framework provides a qualitative understanding of the 
TV^-dependence of infrared-dominated order parameters. In particular, the 
observed Zweig rule violation in the scalar channel could be explained by a 
large variation of the quark condensate from 2_to 3 flavours, close to a chiral 
phase transition. If forthcoming experimentsu related to ir — it scattering 
should soon pin down 2(2), an extensive study of SB^S in the (Nf, N c ) plane 
remains necessary, and could be undertaken through Leutwyler-Smilga sum 
rules joined with lattice computations. 
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